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Abstract. 

We consider classical and quantum propagators for two different time intervals. 
If these propagators follow one another in a Fibonacci sequence we get a discrete 
quasiperiodic system. A theorem due to Nielsen provides a novel conserved quantity 
for this system. The Nielsen quantity controls the transition between commutative 
and non-commutative propagation in time. The quasiperiodically kicked oscillator 
moreover is dominated by quasiperiodic analogues of the Floquet theorem. 

1 Introduction. 

For the general description of time-dependent systems in quantum mechanics we 
refer to pE], [IH], [8], [9], [3]. We consider the time evolution for a particular class 
of time-dependent classical/quantum systems which are periodic or quasiperiodic in 
time. 

We restrict our attention to time evolutions, built piecewise from strings with clas- 
sical symplectic propagators g(t) in 2D phase space. These describe oscillators and 
dilatations, or positive and negative 5-kicks. Quantum propagators are constructed 
as unitary representations S g ^) of the classical linear canonical transformations for 
g{t) from the symplectic group Sp(2,R). 

We present Nielsen's theorem on automorphims of the free group F 2 with two gen- 
erators and from it derive a conserved quantity in terms of the group commutator. 
We analyze Floquet properties of classical and quantum propagators. We show the 
presence/absence of Floquet quantum numbers for periodic system under the control 
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of system parameters. For the periodically kicked classical oscillator we demonstrate 
the occurrence of Floquet bands and Floquet gaps. 

For periodic systems built from two basic intervals, we construct Nielsen's conserved 
quantity as a function of the system parameters. We extend the analysis to quasiperi- 
odic strings. By building the strings through a Fibonacci automorphism we are able 
to survey their quasi-Floquet properties. 



2 Floquet theorem for linear differential equations. 

1 Theorem Floquet: Consider the real first-order periodic linear system of deq. 

— Y(t) = A(t)Y(t), A{t) real, A(t + T) = A(t), T > 0. (1) 

(Jjv 

A fundamental system of two solutions can be written in the form 

Y{t) = exp(±zgt)F (t), Y {t + T) = Y (t), < Q < ^. (2) 

We denote the real number Q as the Floquet index and exp(iQT) as the Floquet fac- 
tor. It will become clear in what follows that the Floquet theorem requires complex 
solutions of the deq eq{T]for real A(t). 

Proof: The discrete time translation T : A(t) — > A(t + T) is a symmetry operation 
of A(t). It generates the abelian group with elements ±nT,n = 0,1,.... The 
unitary representations of this group have the form 

2n 

T -> exp(zfiT), D n (nT) = exp(mfiT), < Q < —. (3) 

The orthogonality and completeness relation of are given in Appendix A. 
Consider a solution of eq. [2] which can be arranged so that it transforms according 
to an irreducible representations of C^, 

Y n {t) :Y n {t + T) = Y n (t)exp(tQT). (4) 

iFrom this solution define the function 

Y^{t) := exp(-^)y n (t) (5) 

Then it easily follows that Y^(t) is periodic, 

Y n (t + T) = Y Q (t) (6) 

Since A(t) is real, the complex conjugate solution Y n := Y~ n is a second linearly 
independent solution with the complex conjugate Floquet factor. With Q = Q we 
obtain eq. [21 
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Remark (i):The Floquet theorem is mathematically equivalent to the Bloch theorem 
for systems with discrete position symmetry, which has the same abstract symmetry 
group Coo. 

Remark (ii): We remark that the periodic system eq. [JJ for certain ranges of its 
parameters cannot have solutions with the Floquet property. We shall give examples 
of this situation in later sections. 

3 The free group F<i and its automorphisms. 

We refer to [TH] for what follows. We also quote [?] as a useful reference to basic 
concepts for free monoids. 

2 Def. The infinite free group F 2 consists of all words generated by concatenation 
from two non-commuting invertible generators {2/1, 2/2}- Its elements can be viewed 

tree formed by words of increasing length. 

3 Def: An automorphism of F 2 is an invertible homomorphic map 

: F 2 -> F 2 , {2/1,2/2} -> {^>i(2/i, 2/2), 02(2/1, 2/2)}, {01,02} e F 2 , (7) 

and so preserves group multiplication and unit element. The set of all automor- 
phisms of F 2 forms the infinite group {0} = Aut(F 2 ). This group, as was shown by 
Nielsen [18], is again finitely generated, with explicit generators given in [TH]. Simi- 
lar results hold true for the free group F n and its group of automorphims Aut(F n ). 
Homomorphic properties of finitely generated groups have the virtue that they can 
be verified by checking them only for the generators. 
Example 1: The Fibonacci automorphism is defined by 

(pfibo ■ {yi, 2/2} -> {<f>fibo,i(yi, 2/2), <f>fibofl(yi, 2/2)} := {2/2, 2/12/2}, (8) 

4>~fKo ■ {2/i>2/2} -> {2/22/r 1 ,2/i}, 

We shall use in particular iterated automorphims as (0/j6 O ) n . Then the words 
(0/ifeo)"(2/i), with their length |iu|, defined [7] as the number of letters in w after 
reducing w by use of all relations of the type yy~ x = e, become 

(<j)fibo) n (yi) ■ 2/1 2/2 2/12/2 2/22/12/2 2/12/22/12/2 ••• /q\ 
KhiboTiyi)] : 1123 5 ... [y> 

The word length K^fiboY (yi)\ is the Fibonacci number rij 

Tij : n j+2 = rij + n j+x , n x = n 2 = 1. (10) 

Notice the recursive structure of the images under ((pfibo) n - Define (0/i& o ) ( ' n ~ 1 - ) (2/i) =: 
Wi, (0/ibo) (n_1) (2/2) =: w 2 . Then it follows from eq. that (0/ih o ) n (2/i) = w 2 , 
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(4>fibo) n {y2) = (^1^2)) an d so we must in step n simply concatenate the two words 
(wi, w 2 ) obtained in step {n — 1). 

4 Nielsen theorem on Aut(F 2 ): Consider the commutator K, = y\y2V\ X V2 V in F 2 
and a fixed automorphism <p(F 2 ) G Aut(F 2 ). Then the image of the commutator 
under <fi obeys, [TH] theorem 3.9 p. 165, 

<j>{K) = wJC^w- 1 , w G F 2 , (11) 

where the sign ±1 and the element w depend on the chosen automorphism (f). 
Proof: It suffices to prove the result for the finite set of generators of Aut(F 2 ) given 
by Nielsen [IE]. 

So the commutator JC is conserved, up to inversion and conjugation, under any 
automorphism of Aut(F 2 ). No counterpart of the Nielsen theorem 3 is known for 
Aut(F n ),n > 3. 

Example 2: Evaluation of the image of the commutator under the particular Fi- 
bonacci automorphism 4>fibo yields 

JC := ymy^yZ 1 , (pfibo(IC) = y 2 (yiy2)y 2 1 (ymy 1 = y2Viy 2 1 yi 1 = /CT 1 . (12) 

The power of the Nielsen theorem appears through the notion of automorphisms 
induced on some group G. 

5 Def: Induced automorphism: Consider a homomorphism horn : F 2 — > G defined 
in terms of the generators by horn, : {yi,y 2 } — > {pi,^} G G. The automorphism 
of G induced by the automorphism (f> G Aut(F 2 ) eq. El is the map 4>g '■ {91,92} — * 
{<pi(gi, g 2 ),<j) 2 (gi, g 2 )}. Clearly induced automorphisms on G are homomorphic to 
the corresponding automorphisms from Aut(F 2 ). 

For an induced automorphism, the commutator induced in the group G has additional 
significance since it measures the non-commutativity of the group. In case JC = I 
we have commutativity on the induced level which allows to rearrange the order of 
successive propagators. 

Example 3: We shall consider in particular automorphisms induced on the symplec- 
tic group Sp(2,R), the group of linear canonical transformations of a ID classical 
hamiltonian system. For g G Sp(2,R), we have the algebraic equivalence relation 
g^ 1 ~ g and the trace relation tr((7~ 1 ) = tr(g). Note that the character is a class 
function. If these properties are combined with the Nielsen theorem, one finds 

6 Prop: For any automorphism of the free group F 2 induced on Sp(2, R), the half- 
trace of the commutator, 

±tr(fc) = itrfotfatfV 1 ) (13) 
is a quantity conserved under iterated automorphisms. 

The condition |x(^) = 1 is a necessary but not a sufficient condition for JC — e. So 
commutativity requires extra checking. 
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If by iterated automorphisms we can generate a string of arbitrary length, the con- 
served quantity eq. [T3] can be assigned to the infinite string. 

If we choose the Fibonacci automorphism eq. Eland associate two matrix propagators 
to its two intervals, we generate a quasiperiodic system. We term such a system 
discrete quasiperiodic. Such system form a special class among general quasiperiodic 
systems [Tfi"] . 

In later sections, we shall associate to the classical system a quantum system. To a 
symplectic matrix g we shall associate a representation by a unitary operator U(g) 
in Hilbert space. 

Since this association is a homomorphic map from symplectic matrices into propa- 
gators which represent linear canonical transformations, the Nielsen theorem then 
extends to a theorem in Hilbert space. 



4 Geometry and dynamics of systems of traces. 



The traces of the symplectic matrices g G Sp(2, R) obviously play an important part 
in the classical and quantum analysis. For a subclass of systems, the computation of 
these traces does not even require the computation of the underlying matrices [I 
The Fibonacci system Example 1 belongs to this class. The trace invariant eq. 
geometrically becomes a cubic surface in 3-space whose shape depends on the value of 
the invariant. The traces under repeated automorphisms and with increasing string 
length determine recursively a discrete dynamical system on this conserved surface. 
We sketch this approach: For the Fibonacci system eq. El we define 2/3 := y\y 2 and 
93 := 9i92- With the induced matrix automorphism written as 



(X n , Y n , Z T 
and the matrix decomposition 



g e Sp{2, R): g = S{g) + A(g), S = -trace(^)e, A = g - S, trace(A) = 



(14) 



(15) 



one finds from eqs. (57-63) of [H] the recursion relations and the conserved quantity 



" S(X(«+!)) " 













rn"\ 



S(K) 



S(Y 
S{Z 

2S{Y n )S{Z n ) - S{X 
4S( gi )S(g 2 )S(g 3 ) + 2 [(S( 9l )S{ gi ) + S(g 2 )S(g 2 ) + S(g 3 )S(g 3 )} - e 



(16) 



All the quantities in eq. [TH] by eq. can be expressed in terms of real traces. The 
recursion relation in eq. HH] determines a discrete dynamical system with initial data 
given by the traces of {gi,g%,gz)- By use of eq. [151 the conserved quantity |x(^) 
can be expressed in terms of the three initial traces. We refer to [IT] for a view of 
the cubic surface which represents ~\;(/C). 
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5 Quadratic hamiltonians and propagators. 



The quantum propagators for quadratic hamiltonians are given in [TP] . [TT] . They 
are unitary integral operators S^rta-ti) where g(t) is a linear canonical transformation 
g G Sp(2m, R), which can be determined by exponentiation from the action of the 
hamiltonian on the canonical position and momentum operators. So for quadratic 
hamiltonians we have a very clear correspondence between classical and quantum 
propagation: The classical time evolution is given by the action of g(t) on the initial 
canonical coordinates at time zero. The quantum propagator is given by the action 
of SgU) on the initial state at time zero. 



6 Hamiltonians and propagators for n = 1. 

In [10], [H] one finds for quadratic hamiltonians the matrix analysis of propagators 
in various dimensions m. The classification of types of hamiltonians for n — 1, 2 is 
discussed in detail in [2T]. We briefly summarize the well-known results for m — 1. 
For m — 1, that is for 1-dimensional position space, we have g(t) e Sp(2,R). This 
group is particularly transparent. Its classes admit a subdivision into three types, 
depending on the value of their half-trace \x(9) := ^ r G?)- Since det(g) = 1, the 
eigenvalues of g are completely determined by the trace. Moreover the trace is inde- 
pendent of symplectic similarity transformations, xiQOQ" 1 ) = x(fiO, an d so allows to 
determine the class of the symplectic matrix g. For I : ||x(flOI < 1, (? is symplectically 
equivalent to the dynamics generated by an oscillator hamiltonian, 

/ : H = -^P 2 + ™X\ (17) 
2m 2 y J 

for // : ||x(sO| > 1 they are equivalent to an inverted oscillator or a dilatation 

II:H = kI(PX + XP) = k(XP+^), (18) 

and for /// : ||x(<?)| = 1 they are, apart from the extra case g = e, |x( e ) = 1, 
equivalent to a free motion, 

III : H = —P 2 , (19) 

2m ' v ' 

All three hamiltonians are hermitian and time- independent. Hence their propagators 
S g (t) are unitary operators. For type I and III the propagators can be written 
as integral operators according to Moshinsky and Quesne [20] in the Schroedinger 
representation or in the Bargmann space of analytic functions, compare jlUj . 
Type I: Harmonic oscillator. 
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Rewrite the harmonic oscillator hamiltonian eq. [T7] as 



H 



1 



-P 2 + mcuX 2 



muo 



(20) 



The corresponding linear canonical transformation is obtained as 

1 











y/mcu 

cos(ut) 
— (mu) sin(W) cos(o;t) 



cos(a;t) sin(a;t) 
— sin (cut) cos (cut) 

— sin(a;t) 



^Jmuj 

~r 



(21) 



Eq. [21] admits the complex diagonalization 



cos(a;t) sin(a;t) 
— sin(a;t) cos(u;i) 



(22) 



1 i 
1 -i 



exp(—iut) 

exp(iuit) 



1 

—i 



In the complex domain we can therefore pass to a pair of classical variables which 
propagate with exp(^fiut). 

The corresponding quantum propagator in the Bargmann space of analytic functions 
is [ID] 

(z'\Sg c (t)\z) = exp(— iut/2) exp(exp(iut)z'z). (23) 

The propagator in the Schroedinger representation is given in [TU] pp. 34-6. 
Type IT. The hyperbolic and dilatation hamiltonian. 
For type II the symplectic matrix may be taken as 



g(t) 



exp(/tt) 

exp(— nt) 



(24) 



Note the correspondence to eq. [221 for K — > — iuo. 

There is for type II no obvious integral operator representation. In this case we 
proceed directly from the hamiltonian eq. [18] and find 



i t d 

exp(—Kt—H)ip(x) = exp(— «-)(exp(— ntx— )ip){x) 

11/ £ (JjJb 



exp(— K-)ip(exp(— nt)x) 



(25) 

This propagator is a time-dependent unitary dilatation. Applied to the ground state 
of an oscillator with frequency c^ch it transforms it into an oscillator ground state with 
frequency u(t) = u exp(—2Kt). As a function of time, the uncertainty in position 
would increase while the uncertainty of momentum would decrease. 
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Another form for class // is given by the hyperbolic hamiltonian 



H = - 
2 



p2 

mK 2 X 2 

m 



(26) 



with the symplectic propagator 
9«(t) = 



cosh(ftt) (m/c) 1 sinh(/?t) 

(mK) sinh(Kt) cosh(Kt) 



(27) 



By choosing at time t = the initial matrix 



h(0) 



1 

rriK 



det(h(0)) = mK, 



{21 



we determine two canonical column solutions given by 

= g K (t)h(0) = 



X z (t) X n (t) 
Pi(t) P n (t) 



cosh(/tt) sinh(/tt) 
(mK) smh(Kt) (mK) cosh(/tt) 



(29) 



Inserting the fundamental two solutions into the hamiltonian eq. [2S we find the two 
different energies 

2 2 
„ TTIK „ ,„v TfLK . > 

Ei(0) = - — , E n (0) = —. (30) 
The types I, II of propagation in phase space are shown in Fig. 1. 





Fig. 1. Phase space representation. The harmonic and the hyperbolic column 
solutions I: eq. EH and II: eq. [27] in phase space. 

In appendix B we give the relation between the propagator and the transfer matrix. 
Type III: Free particle. ^From the hamiltonian 

H = —P 2 (31) 
2m v ; 
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one finds [TU] the symplectic matrix 



g(t) 



1 -t 
o 1 



(32) 



and the propagator 



(x'\S g ^\x) = exp(-i7r/4) 



2nht 



m 



exp(i 




(33) 



7 Prop: The types I, II, III together with the case g = e exhaust, up to conjugation 
within Sp(2,R), the symplectic propagators with quadratic hamiltonians. 

7 Quadratic piecewise time-independent hamilto- 
nians. 

We consider now a system whose hamiltonian for < t < T = T\ + T 2 is time- 
independent on the two intervals. Assume that the propagation for the full interval 
T is repeated periodically. Over the period we have then a propagation with the 
half-trace given by eq. H2J 

For < 1 we can write \x = cos(fiT) and determine Q, < Q < 2ix/T. This value 
Q is the Floquet index of a classical symplectic periodic system. The two Floquet 
factors exp(±zf2T) appear as the eigenvalues of the symplectic matrix gi(Ti)g 2 (72) 
and determine the classical propagators U(T). The Floquet index is completely 
determined by the parameters (u\, T\, u; 2 , T 2 ) of the system. At multiples nT,n = 
1,2,... the classical propagator becomes (U(T)) n . 

In the quantum propagation, Q represents the overall frequency of an effective oscil- 
lator hamiltonian. The propagator over the period T from eq. [23] becomes 



We emphasize that this result holds true only for T = Ti+T 2 whereas at intermediate 
time the propagator is governed by giif) or g 2 (t) respectively. The classical Floquet 
parameter enters the quantum propagator in an essential way. The eigenstates \N) 
of the overall harmonic oscillator of frequency Q from eq. [3H are the monomials of 
degree N. These eigenstates under the discrete time translations transform according 
to 



This again generates an irreducible representation of the time translation group. By 
the unique decomposition NQ = m27i/T + Q',0 < Q' < 2ir/T,m = 0,1,2..., this 
representation becomes D (NT) = D n (T). 



(z'\Sg c (T)\z) = exp(— iflT/2) exp(exp(iQT)z' z), 



(34) 



T : \N) -»• \N)exp(iNQT) 



(35) 
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It follows that not the quantum propagator itself, but only the eigenstates in its 
decomposition transform irreducibly under the period T. 

For ^\x\ > 1 the situation changes. By writing ||%| = cosh(/tT) we determine an 
effective classical dilatation parameter k. The eigenvalues exp(±KT) represent this 
dilatation over one period. The corresponding classical propagation does not repre- 
sent a Floquet propagation since the eigenvalues of the symplectic matrix do not have 
absolute value 1. They form a non- unitary real representation of the time translation 
group. The nth powers of these two factors now approach 0, oo respectively. 
The unitary quantum propagator over a period T from eq. [25] becomes 

% T d T 

exp(—K,T—H)ip(x) = exp(—K—)(exp(—K,Tx—)i(j)(x) = exp(— k— )ip(exp(— k,T)x). 
it 2 dec 2 

(36) 

Again, this propagator is valid only at time T = T\ + T 2 . For intermediate time the 
propagator is governed by one of the oscillator matrices gi(t) or g 2 (t) respectively. 
The quantum propagator over the periods T stays unitary, but for n — > oo when 
acting on a state transforms it up to normalization into a 5-distribution of sharp 
effective momentum/position. 

Remark (i): The results on time-periodic systems from piecewise quadratic hamil- 
tonians resemble the analysis of position-periodic Kronig-Penney systems p2] with 
piecewise constant potentials. In the latter systems, the potentials play the role of 
the frequencies, and the Bloch /c-label plays the role of the Floquet index. Real 
values of the Bloch parameter yield a band index, while imaginary values indicate 
exponentially increasing or decreasing solutions of the Schroedinger equation which 
are excluded as band gaps. 

Example 5: For the system of section 8 in case \ \x\ < 1 we obtain an overall Floquet 
state with Q determined by \\ = cos(fiT). Each Floquet state is the analogue to a 
Bloch state of a band model in x-space. 



8 Oscillators and 5-kicks. 

The transfer and propagator matrix methods apply to sequences of time intervals. 
Example 6: The classical square well for the time interval t has the type / oscillator 
propagator 

cos(co't) (mw)" 1 sin(o;t) 

— (moo) sin(fca) cos(a>t) 



9u,(t) 



(37) 



The initial values for the two fundamental column solutions /, II we take as 



Mo) 



1 

muj 



(38) 
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Then the energies computed by insertion into the harmonic oscillator hamiltonian 
become 



Ej = E 



11 



(39) 



Example 7: The square tunnel for the time interval t has the repulsive oscillator 
type II propagator eq. [27J The energy of the two fundamental solutions are Ej = 



ran" rn ran' 

— ,E n - — 



Following what is done in space in [T3] , [Ej we introduce w.r.t. time, negative and 
positive £-kicks. They are obtained as limits of square well and square tunnel intervals 
of finite length. In the limits : muj 2 T — > u', ujT — > and mn 2 T —>■ u', kT — > with 
u the strength parameter we get from eqs. [37J, EH] respectively the symplectic transfer 
matrices for negative/positive 5-kicks in the form 



g±u' 



1 

Tu' 1 



(40) 



9 The periodically kicked oscillator. 



We employ here the transfer matrix, see eq. 
M(t) = Xi = 



with u 1 



cos(a;Tj) uj 1 sin(a;Tj' 
— uj sin(a;Tj) cos(a;Tj) 



= mti, 

1 " 

u 1 



1,2. 



(41) 



with positive 5-kicks. The half-trace becomes 



1 u 
X = -Tr(M(T)) = cos(^T) + — sin(^T) 



This expression has the period — . We find 

X = : cot (u;T ) 



u 



(42) 



(43) 



The solutions To of this equation exist in each interval 2mir < uoT < 2(m+ 1)tt, m = 
0, ±1, ±2, ... and mark the centers of Floquet bands as functions of k. The equations 



^X = +1 : W^±; 



•1 : cot( 



u 

u 



(44) 



determine for any such band its edges T± as functions of k. Note that these equations 
have in ujT± the period Atc. In between two Floquet bands there are gaps. In these 
gaps we have ||%| > 1 and so solutions of type II. 
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7 Prop: The Floquet spectrum of the periodically kicked oscillator. The 

periodically kicked oscillator has a Floquet band spectrum. Each interval Im-n < 
ujT < 2(m + 1)tt, m — 0, 1, 2, .. carries two bands. The centers To from eq. 03] are 
repeated periodically, the edges given by eq. HU The bands are separated by gaps 
which admit only solutions of type II. 

Example 8: A simple example is provided by the choice ^- = 1. Then the half-trace 
\x ec L- H2]with := toT becomes 

I x = cos(/3) + sin(/3) = \/2cos(/3 - vr/4). (45) 

The band edges + ,0~ and the band centers 0° are located at 

+ : + 27rmi, | + 27rm 2 , (46) 

37T 

: 7r + 27rm 3 , — + 27rm 4 , 

5° : - + 7rm 5 , m,- = 0, ±1, ±2, ... 
4 

Bands and gaps alternate with the same width |, as shown in Fig 2. 



+ 1 



-1 




CUT 

Fig. 2. Floquet band structure for the periodically kicked oscillator. Shown 
is the half-trace \x eq. H2] as a function of ujT for the parameter j*- = 1 of eq. HS1 
Bands and gaps alternate periodically and both have the width A0 = |. 



12 



10 The quasiperiodically kicked oscillator. 

For properties of systems with a general quasiperiodic Hamiltonian we refer to [5]. 
The Nielsen theorem is valid for all discrete quasiperiodic systems of Fibonacci type, 
but its significance can already be seen on a simple example. Here we study a specific 
discrete quasiperiodic system of Fibonacci type in time, in analogy to a system in 
space from [12J. Following |6J, it would be worthwile to study the spectral function 
for this discrete quasiperiodic system. 

For two intervals (1, 2) which form a quasiperiodic Fibonacci sequence we take the 
matrix propagators in the form 



A; 



cos(o;Ti) uj 1 sin(a;Tj) 
— uj sin(a; Tj) cos(a>Tj) 





' 1 


" 




u 


1 



1,2. 



(47) 



Then we compute, expanded in powers of ^, the products 
A1A2 = 



(4 



u 
+ - 

UJ 



+( 



.It. 



uj 



cos(a;(Ti + T 2 )) uj 1 sin(a;(T 1 + T 2 
-uj sin(co>(Ti + T 2 )) cos(u;(Ti + T 2 )) 

2 sin(a;T 1 ) cos(a;T 2 ) + cos(u/Z~i) sin(a;T 2 ) uj' 1 sin(a;T 1 ) sin(u;T ; 
uj{2 cos(a;Ti) cos(a;T 2 ) — sin(u/Z~i) sin(a;T 2 )) cos(a;Ti) sin(a;T 2 ) 

sin(a;Ti) sin(u;T 2 ) 
uj cos(u;Ti) sin(c<jT 2 ) 

cos(w(Ti + T 2 )) -uj' 1 sin(cj(Ti + T 2 
ujmi(uj(Ti + T 2 )) cos(co>(Ti + T 2 )) 

sin(a;Ti) cos(a;T 2 ) — uj' 1 sin(a;Ti) sm(ujT 2 ) 

uj(2 cos(a;Ti) cos(a;T 2 ) — sin(cjTi) sin(cjT 2 )) 2 cos(u/Z~i) sin(a;T 2 ) + sm(ujTi) cos(u/T 2 ) 


—uj sin(u/Zi) cos(a;T 2 ) sin(a;T 1 ) sin(o;T 2 

For the commutator JC and its inverse we find 



(A 2 A: 



u 



+- 




UJ 








M- 




UJ 





/C = (A 1 A 2 )(A 2 A 1 )- 1 



(49) 



1 
1 



u 



0J 



sin(a)(Ti - T 2 )) 



7/ 

+ (-) 2 sin(^(T 1 -T 2 )) 

UJ 



cos(a;(Ti + T 2 )) -uj 1 sin(co>(Ti + T 2 )) 
-uj sin(tu(T 1 + T 2 )) - cos(u;(T 1 + T 2 )) 

sin(a;Ti) cos(a;T 2 ) —uj' 1 sin(k/Z~i) sin(u/T 2 ) 
uj cos(a;T 1 ) cos (u;T 2 ) — cos(a;T 1 ) sin (u;T 2 ) 



KT 



(A 2 A 1 )(A 1 A 2 )- 1 



1 
1 



u 



+ -sm(uj(T 1 -T 2 )) 



UJ 



+ ( 



. u , 



UJ 



sin(a;(T 1 - T 2 )) 



- cos(u;(Ti + T 2 )) uj 1 sin(cj(Ti + T 2 )) 
u; sin(u;(Ti + T 2 )) cos(tu(Ti + T 2 )) 

cos(u/Z~i) sin(o;T 2 ) uj' 1 sin(u;Tx) sin(u;T 2 ) 
uj cos{ujTi) cos(u;T 2 ) sin(u/Z~i) cos (u/T 2 ) 
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Remember that in each iteration of the Fibonacci automorphism the commutator 
/C by eq. [T2] is transformed into its inverse. This does not affect the trace which 
becomes an invariant under the automorphism. 
The half-trace of the commutator therefore becomes the invariant 

X = ~Tr(X:) = -Tr^" 1 ) = 1 + ^HsinMT! - T 2 )) 2 . (50) 

Of interest are the points where the invariant takes the value 7 = 1. These are given 
from eq. E01 by 

mi(u(T 1 -T 2 )) = 0, (51) 
(jj(T\ — T 2 ) = /mi", m = 0, ±1, ±2, .... 

For the Fibonacci case T\ = T,T 2 = tT we find from eqs. [501 EH 

cu(t-1)T = mw,m = 0,±1, ±2, .... (52) 

The full commutator K, at these points from eq. H9] becomes 

K. = I. (53) 

and we find 
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8 Theorem: The invariant I as a function of T 1; T 2 has the value X = 1 at the points 
eqs. EH and E21 These points are repeated periodically. At each of these points, the 
commutator K, becomes the unit matrix and so the two matrices eq. H3 commute. 

For the Fibonacci model of eq. E21 we get periodic points where the two matrices 
Ai,A2 commute. This does not automatically imply that we have a quasiperiodic 
Floquet propagation. To assure that, we must for both matrices be in an individual 
Floquet band. The relevant parameters for the individual Floquet bands are uT, tujT 
respectively which scale by r. To find overall propagation we must superimpose the 
two band structures and find the regions of band overlap. Moreover for commuta- 
tivity we must fulfill the condition eq. E2J The superposition of the bands and the 
location of commutativity are shown in Fig 3. 






5 

t 


10 


15 

t 


20 


25 


30 

t 


WMMVAMM/ 



Fig. 3. Quasiperiodic Fibonacci system. Bottom: For the parameters of eq. 
351 the half-traces of the two systems exhibit two systems of bands as functions of 
uoT. Top: Overlap regions of the two systems of bands and the invariant X{ujT). 
If a commutative point with X = 1 hits a band overlap region, the system shows 
quasiperiodic Floquet behaviour. These points are marked by vertical arrows. 

The propagation in phase space for multiple values of a commutative in-band case 
is shown in Fig. 4. All the points reached fall on a single ellipse equivalent to a 
harmonic oscillator. 
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Fig. 4. Quasiperiodic Floquet theorem: For multiples of a value uT with X{uoT) = 1, 
hence commutative, and moreover in the overlap of the two Floquet bands, all phase 
space propagators fall on an ellipse equivalent to a harmonic oscillator. 

11 Conclusion. 

The algebraic Nielsen theorem applied to classical and quantum quasiperiodic sys- 
tems yields a novel conserved quantity. This quantity as a continuous function of 
system parameters controls transitions between commutative and non-commutative 
propagation in time. The propagation of the quasiperiodically kicked oscillator 
demonstrates an interplay between the Nielsen and Floquet properties. 
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Appendix A: Irreps of the discrete translation group. 



We give here the orthogonality and completeness relations for the discrete translation 
group Coo as discussed in [15]. In the present case they apply to Floquet states. 
The orthogonality relation is 



D n (nT) = exp(27riQnT), 

Vie BZ = {0, 2n/T}, n e {0, ±1, ±2, ...} 

^lP(nT)D n ' (nT)dfl = \BZ\5{Q - Of). 



(54) 



Eq. [5H yields an orthogonality relation between evolutions belonging to two different 
Floquet indices. It does not imply an orthogonality between pairs with the same 
Floquet index from different Floquet bands, as they are found for example in section 
10. The completeness relation is 



neBZ 



D n (nT)D Q (mT)dn = \BZ\5 n , r 



(55) 



Appendix B: Classical propagator and transfer ma- 
trix. 

The matrix which propagates the canonical pair of position and momentum {X(t), P(t)} 
we denote as the classical (symplectic) propagator g(t) e Sp(2,R), det(g{t)) = 1. 
The matrix which propagates a fundamental pair of position and velocity {X(t),X(t)} 
we denote as the transfer matrix M(t). Again we may choose det(M(t)) = 1. For a 
hamiltonian with standard kinetic energy we have the relation P(t) = mX(t). ^From 
this we get the conjugation relation 



M(t) 



*m 

o 4= 



g(t) 



-7= 




m 



(56) 



between the symplectic propagator g(t) and the transfer matrix M(t). From the rela- 
tion eq. [56] it follows that for a sequence of strings both their symplectic propagators 
and their transfer matrices are multiplied in the same order. 

The first and second column of the propagator matrix g(t) each yield a fundamental 
canonical pair {/, II} of solutions and trajectories for the hamiltonian equations of 
motion. For quadratic hamiltonians, the hamiltonian equations of motion are linear. 
Then the most general classical solution is a linear superposition of fundamental 
solutions. The superposition changes the classical energy which therefore must be 
recalculated for any superposition. 
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